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Abstract. We report the discovery of an eclipsing binary 
- HS 0705+6700 - being an sdB star with a faint compan- 
ion. From its hght curve the orbital period of 8263.87s, 
the mass ratio of the system q — 0.28, the inclination of 
84?4 and other system parameters are derived. The com- 
panion does not contribute to the optical light of the sys- 
tem except through a strong reflection effect. The semi- 
amplitude of the radial velocity curve Ki = 85.8 km s~-^ 
and a mass function of /(m) = 0.00626 are deter- 
mined. A spectroscopic analysis of the blue spectra re- 
,—1 [ suits in Toff = 28 800 K, log g = 5.40, and log(nHc/"^H) = 

■ —2.68. These characteristics are typical for sdB stars, as is 
' its mass of 0.48 Mq. According to its mass (0.13 Mq) and 
• radius (0.19 i?©), the companion is an M dwarf. The pri- 

mary is in a core helium burning phase of evolution, and 

■ the system must have gone through a common envelope 
I stage when the primary was near the tip of the red giant 

branch. 

^ ■ Key words: subdwarfs - binaries: eclipsing - binaries: 
I spectroscopic - stars: early-type - stars: fundamental pa- 
rameters - stars: individual: HS 0705+6700 
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1. Introduction 

Subluminous B stars (sdB) dominate the populations of 
faint blue stars of our own Galaxy and are found in both 
the old disk and in halo populations, e.g. as stars form- 
ing the blue tails to the horizontal branches of globular 
clusters (Ferraro et al. 1997). 

In the context of galaxy evolution sdB stars are impor- 
tant because they are sufficiently common to be the domi- 
nant source for the "UV upturn phenomenon" observed in 



elliptical galaxies and galaxy bulges (Brown et al. 1997, 
20004 Greggio & Renzini |1990| , |1999D . They might also 
be used as age indicators for elliptical galaxies (Brown et 
al. |2000b|). 



However, important questions remain over the evolu- 
tionary paths and the appropriate timescales. There is 
general consensus that the sdB stars can be identified 
with m odels for Extreme Horizo ntal Branch (EHB) stars 
(Heber |l986i Saffer et al. |l994|). Like aU HB stars they 



are core helium burning objects. However, their internal 
structure differs from typical HB stars, because their hy- 
drogen envelope is very thin (<1% by mass) and therefore 
inert. As a consequence EHB stars evolve directly to the 
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white dwarf cooling sequence thus avoiding a second red 
giant phase (Dorman et al. 1993). 

The discovery of multi-mode, short-period (P 2 — 10 
min.) pulsators among sdB stars (see O'Donoghue et 



al. 1999 for a review) has opened a new attracting pos- 
sibility of probing their interiors using seismological tools. 
Recently, Brassard et al. (2001) were able to derive the 
mass and hydrogen envelope mass of the pulsating sdB 
star PG 0014-1-067 by asteroseismology. Their results 
(total mass M = 0.49 ± 0.019 and envelope mass 
log(A/onv/ ^^o) = —4.3 ±0.22) are in excellent agreement 
with the predictions by evolution theory. 

Considerable evidence is accumulating that many sdB 



stars reside in close binaries (Maxted et al. 2001a, Saffer et 
al. 2001). Therefore mass transfer should play an impor- 
tant role in the evolution of such binary systems. Detailed 
investigations of sdB binaries, in particular eclipsing sys- 
tems, are crucial to determine their masses. However, only 
two such eclipsing binaries, HW Vir (Menzies & Marang 



1986) and PG 1336-018 (Kilkenny et al. 1998) are known 



up to now, which consist of an sdB star and an optically 
invisible M dwarf companion. Here we report the discov- 
ery of the third such system ~ HS 0705-^6700. 

HS 0705-1-6700 was selected from a list of candi- 
dates drawn from the Hamburg Schmidt survey (Hagen 



et al. 1995). Follow-up spectroscopy (Heber et al. 1999 
Edelmann et al. 20011j ) revealed that its effective tem- 



perature lies in the domain predicted for the pulsational 
instability. Therefore, HS 0705-1-6700 was included in a 
photometric monitoring programme at the Nordic Opti- 



cal Telescope (see 0stensen et al. 2001a, 2001b) in order 



to search for pulsations. Indeed, the star turned out to be 
variable, but not on the time and flux scale expected for 
pulsations. The onset of a primary eclipse event shortly 
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after the start of the observations brought the observed 
brightness down by about a magnitude. This made it evi- 
dent that HS 0705+6700 must be an echpsing binary. This 
was confirmed in the foUowing night when both primary 
and secondary minima were monitored. 



2. Observations 

2.1. Photometry 

The initial observations were made on 5 and 6-Oct-2000 
at the Nordic Optical Telescope (NOT) with the Andalu- 
cia Faint Object Spectrograph and Camera (ALFOSC), 
equipped with a Loral, Lesser thinned, 2048x2048 CCD 
chip, and modified with our own control software to be 
able to observe in high-speed multi- windowing mode. The 
sky area available for locating a reference star is limited 
to the area of the chip: ~ 6.5x6.5 arcminutes^. On 5-Oct- 
2000, the observations started at UT 05:38 and covered 
a time period of 3590 s, and on 6-Oct-2000 measurements 
started at UT 03:28 and extended over 10280 s. 

We observed HS 0705+6700 using three reference stars 
for constructing the relative light curve. The observations 
were made with a Bessel S-band filter in order to optimize 
for detection of low level pulsations. A cycle time of 20 s 
was used, allowing an actual integration time of 16.5 s. A 
total of 634 exposures were collected in the nights of 5 and 
6-Oct-2000. The data were reduced online using the Real 
Time Photometry (RTF) program developed by one of us 
(RO) as part of his PhD project (Ostensen |2000D . More 
details about this software are given elsewhere (Ostensen 
et al. ^OOlal ). 

The i?-band photometry of HS 0705+6700 was ob- 
tained on 2-Nov-2000 at the Calar Alto Observatory using 
the 2.2m telescope with the CAFOS instrument. Obser- 
vations started at UT 01:57, and the total length of the 
data set was 12841s. A total of 598 measurements were 
made with a cycle time of 21 s (the integration time was 
10 s). Each exposure, taken through a Johnson R filter, 
is a 2x2 binned subframe of the 2kx2k blue-sensitive, 
back-illuminated SITe chip. The field around the target 
contains several other stars, five of which were used as 
reference stars. The basic data reduction as well as the rel- 
ative aperture photometry was performed using the IDL 

|l99g| ). The 1 a error of the 



software TRIPP (Schuh et al. 
resulting normalized differential light curve, derived from 
the variations of all relative reference light curves, is 0.013. 

Following the B and R light curves obtained in Oc- 
tober and November 2000, additional photometric obser- 
vations were collected with various smaller telescopes (on 
Tenerife, in Norway, Poland, and Greece) in five differ- 
ent nights scattered over the next five months (November 
2000 to March 2001). Nine more primary minima were 
covered. The measurements were made with CCD cam- 
eras in different filters (i?, /), and in one case with a clas- 
sical photo-multiplier tube (PMT) photometer in unfil- 



tered light. Dates and other details of this complementary 
photometry are contained in Tabic ^. 

Table 1. Spectroscopic observations (heliocentric Julian 
Date for mid-exposure, length of exposure, and heliocen- 
tric radial velocity. 



nel.JJJ (^miaj 


exposure 


nel. HV 


-2451900 


time (s) 


(km s^i) 


79.55107 


900 


56.5 


79.56562 


900 


-18.1 


79.57960 


900 


-74.5 


79.59264 


900 


-110.1 


80.46672 


600 


-114.5 


81.31732 


600 


-120.2 


81.32809 


600 


-106.4 


81.34120 


900 


-38.2 


81.35216 


600 


4.6 


81.36139 


600 


52.5 


81.37098 


600 


32.9 


81.38026 


600 


4.9 


81.39023 


600 


-45.7 


81.39943 


600 


-81.6 


81.40900 


600 


-116.6 



2.2. Spectroscopy 

Optical spectra of HS 0705+6700 were obtained in March 
2001 at the DSAZ observatory at the 3.5m telescope 
with the Twin spectrograph in dichroic mode covering the 
wavelength ranges 3925A to 5015A in the blue channel and 
5975A to 7070A in the red channel at a spectral resolution 
of lA. Details of the observations are given in Table |l|. The 
data were reduced with the MIDAS package distributed 
by the European Southern Observatory ESO. 



3. Photometric analysis 

3.1. Ephemeris 

Our B and R CCD photometry of October-November 
2000 provided clear evidence that HS 0705+6700 is a 
short-period eclipsing binary. The complete phase cov- 
erage and high time resolution of the CCD light curves 
allowed for an accurate determination of minimum times 
and a period analysis. 

The CCD B measurements carried out with the 2.5m 
Nordic Optical Telescope covered a complete primary 
eclipse in each of the two nights of October 5 and 6, 2000. 
The CCD R photometry obtained with the CAFOS in- 
strument attached to the Calar Alto 2.2 m telescope on 
November 2, 2000, continuously extended over more than 
one complete orbital cycle and covered two subsequent 
primary minima. Already from these four minimum times 
we derived an orbital period of 8263.88 seconds. 
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Table 2. Times of primary minima derived from CCD and PMT photometric observations 



Date 


rljij 


Error 


Epoch 


u - 


- o 


MetnoQ 


Filter 


iype 


Telescope^ 


Observers'^ 


05-10-00 


2451822.75978 


0.00005 





-0 


00004 


Kwee 


B 


CCD 


NOT 2.5m 


R0, JES 


06-10-00 


2451823.71648 


0.00010 


10 


+0 


00019 


Kwee 


B 


CCD 


NOT 2.5m 


R0, JES 


02-11-00 


2451850.59301 


0.00011 


291 


+0 


00001 


parabola 


R 


CCD 


CA 2.2m 


SS, JD 


02-11-00 


2451850.68866 


0.00012 


292 


+0 


00002 


parabola 


R 


CCD 


CA 2.2m 


SS, JD 


30-11-00 


2451878.71316 


0.00003 


585 


+0 


00005 


Kwee 




PMT 


lAC 0.8m 


FJ 


09-01-01 


2451919.2670 


0.0002 


1009 


-0 


00029 


Kwee 


I 


CCD 


Skb 0.6m 


sz 


09-01-01 


2451919.3627 


0.0010 


1010 


-0 


00024 


parabola 


I 


CCD 


Skb 0.6m 


sz 


09-01-01 


2451919.4586 


0.0001 


1011 


+0 


00002 


Kwee 


I 


CCD 


Skb 0.6m 


sz 


15-02-01 


2451956.4738 


0.0013 


1398 


-0 


00004 


parabola 


R 


CCD 


Crc 0.5m 


sz 


16-02-01 


2451956.5697 


0.0011 


1399 


+0 


00022 


parabola 


R 


CCD 


Crc 0.5m 


sz 


17-02-01 


2451957.5274 


0.0002 


1409 


+0 


00145 


Kwee 


R 


CCD 


Crc 0.5m 


sz 


27-03-01 


2451996.3586 


0.0001 


1815 


+0 


00011 


Kwee 


R 


CCD 


Kry 1.2m 


sz 


28-03-01 


2451997.3150 


0.0001 


1825 


+0 


00004 


Kwee 


R 


CCD 


Kry 1.2m 


sz 



^ Relative to ephemeris (1) 

^ NOT: Nordic 2.5m, La Palma; CA: Calar Alto 2.2m, Spain; lAC: lAC 0.8m, Teide, Tenerife; Skb: Skibotn 0.6m, Norway; 
Crc: Cracow Jagiellonian Obs. 0.5m, Poland; Kry: Kryonerion 1.2m, Greece 

R0: Roy Ostensen; JES: Jan-Erik Solheim; SS: Sonja Schuh; JD: Jochen Deetjen; FJ: Frank Johannessen; SZ: Stanislaw 
Zola 



For a more refined period analysis nine more pri- 
mary minima were used, which had been observed with 
several smaller telescopes between November 2000 and 
March 2001 (see Table ||). Including these minima means 
a considerable broadening of the time base, allowing for 
more accurate period determination. The minimum times 
were d etermined either by the Kwee- van Woerden method 
( 1956 ) or by fitting parabolas to the cores of the minima. 
Table ^ lists the minimum times together with other ob- 
servational details. 

The standard deviation from a linear O — C represen- 
tation amounts to 1.5 • 10"'* days (« 13 sec), which is less 
than the mean time resolution between subsequent CCD 
exposures when monitoring the light curves (mostly about 
20 sec). Eight minima lie within a, 1 a belt, and four oth- 
ers within a 2 a belt. There is only one primary minimum 
(the one of February 17, 2001, at epoch 1409) deviating 
by about 10 a from the linear regression. This timing was 
therefore not used for the period determination. From the 
remaining 12 primary minima distributed over a time span 
of nearly 6 months, the following linear elements for the 
heliocentric primary minimum were derived: 

HJD 2451822"?75982 + 0'?09564665 • E 
±22 ±39 

This ephemeris was applied for the phasing of the light 
curves and of the spectroscopic observations analyzed in 
this paper. 

3.2. Light curve solution 

The B and R light curves obtained with the 2.5m NOT 
and 2.2m Calar Alto telescopes, respectivel y, wer e ana- 
lyzed with the MORO code (Drechsel et al. |l995|) . This 



light curve solution program uses the Wilson-Devinney lo- 



gistical approach (1971), but is based on a modified Roche 
model for the binary structure, taking into account the ra- 
diative interaction between the components of hot, close 
binaries. Also, the differential corrections procedure for 
parameter optimization is replaced by the more powerful 
■simplex algorithm, which was first appli ed to eclipse light 
curve solutions by Kallrath & Linnell ( 1987 ). For more 
details and applic ation e xamples see, e.g., Drechsel (2000) 
or Drechsel et al. ( 1995 ). 

The B and R light curves were represented by normal 
points formed by binning the fluxes of individual measure- 
ments over narrow time intervals. One reason for this was 
to reduce the observational scattering to a smoother run, 
and on the other hand the computational performance was 
improved. The phase resolution of the normal points was 
generally chosen as 0.005 phase units (41 sec), except for 
the core of the primary minimum, where a step width of 
0.002 phase units (16 sec) was used. The total binning 
window size was 0.02 phase units (165 sec). Hence the 634 
and 598 individual measurements of the B and R curves, 
respectively, were reduced to 212 normal points in each 
passband. Typically, 12 individual observations were aver- 
aged to form a single normal point. The input light curves 
(in intensity units) were normalized to unity at quadrature 
phase 0.25. 

We used the principal Wilson-Devinney mode 2, which 
poses no restrictions on the presumed system configura- 
tion (detached or semi-detached), and links the luminos- 
ity of the secondary component to its effective tempera- 
ture by means of the Planck law. The latter aspect is of 
nearly no importance here, because the luminosity ratio 
between primary and secondary flux contributions is so 




Fig. 1. B and R light curves of HS 0705+6700 obtained on Oct. 5/6 {B) and Nov. 2, 2000 (i?); the sohd hnes 
represent the best fit simultaneous solution according to the parameters listed in Table ^; the R light curve is offset 
by 0.1 intensity units relative to the B curve. 



huge in this system (L1/L2 is of order 10** for B and 10"^ 
for R passbands), that the secondary is photometrically 
only evident through the reflection effect and shadowing 
of primary light. The B and R light curves were solved 
simultaneously to yield consistent solution parameters. 

Since the total number of light curve parameters is 
considerably large (12 + 5 • n), where n is the number of 
spectral passbands (filter light curves), i.e. equals 22 for 
our case of a simultaneous B and R fit, it is important 
to reduce the free parameter set by consideration of spec- 
troscopic and theoretical boundary conditions and other 
consistency checks. Otherwise, too many free parameters 
would easily tend to produce underdetermined solutions 
with no guarantee of uniqueness. 

As was obvious from the sinusoidal shape of the ra- 
dial velocity curve (Sect. p| ) and from the position of 
the secondary minimum at exactly phase 0.50, it was im- 
plicitly assumed that we are dealing with circular orbits 
(e = 0) of synchronously rotating components in this very 
close binary system, for which an extremely short syn- 
chronization time scale of a few decades is expected (Zahn 



1977). According to the early spectral type of the primary 

1; also its gravity 
1 as expected for 



we assumed a bolometric albedo Ai 
darkening exponent was fixed at gi 
radiative outer envelopes (von Zeipel 1924 ). For the cool 



convective secondary, 52 was set to 0.32 according to Lucy 
( 1967 ) . Already after the first trial runs with the secondary 
albedo A2 as a free parameter, it became obvious that the 
enormous reflection effect evident as a broad orbital hump 
centered on secondary eclipse, with an amplitude of about 
12% of maximum light, requires a "mirror-like" surface 
of the tidally locked secondary in the heated area facing 
the primary. Reasonable representations were only possi- 
ble for complete reradiation (energy conservation), i.e. A2 
was fixed at 1.0 in all further runs. Linear limb darken- 
ing coefficients xi{B, R) for the primary were interpolated 
from the tables of Wade & Rucinski ( 1985| ) and fixed at 
xi{B) = 0.26 and xi{R) — 0.19. Since it was known from 
an analysis of similar sdOB binaries with strong reflection 
effect (Hilditch et al. 1996) that the limb darkening co- 
efficients of cool stars irradiated by a hot companion can 
strongly deviate from "normal" values of very cool dwarf 
stars, it was decided to treat X2 as an adjustable parame- 
ter. From our spectroscopic analysis (Sect. ^ it was clear 
that the sdB primary must have an effective temperature 
close to 30 000 K. Though Ti was always used as an ad- 
justable quantity, it actually did not vary much during 
numerous optimization runs, and mostly remained within 
a narrow range between about 29 000 and 30 000 K, com- 
patible with the spectroscopic error margin. 
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Fig. 2. Meridional intersection of the equipotential struc- 
ture of HS 0705+6700 corresponding to the best fit pho- 
tometric sohition for a mass ratio q = 0.278; the inner 
critical Roche lobe is shown as dashed line. 

The remaining set of adjustable parameters comprises 
the orbital inclination i, the mass ratio q = M2/M1, the 
Roche potentials fli and fl2 of the two stellar surfaces, 
the color-dependent luminosity Li of the primary and of 
a possible third light contribution ^3, the effective temper- 
atures Ti and T2 , and the radiation pressure parameter Si 



(see Drechsel et al. 1995); S2 was set equal to zero because 
the secondary temperature is low (about 3 000 K), which 
means that radiation pressure forces exerted by this star 
are negligible. Finally, L2 was not adjusted as an inde- 
pendent parameter, but recomputed from the secondary 
radius r2 and its effective temperature. 

A convergent simultaneous solution of the B and R 
light curves was achieved after several trial runs using 
different start parameter sets, most of which lead to an 
essentially identical final solution. Table |^ gives a list of 
the fixed and adjusted parameters. The corresponding fit- 
ted light curves together with observations are shown in 
Fig. |l|. It is obvious from the graph that the overall run 
of the B and R light curves is very well reproduced by 
the best fit theoretical curves (solid lines). Especially the 
shape of the primary eclipse minima and the pronounced 
reflection effect are fitted very well. A minor systematic 
deviation is apparent in the secondary minima: the fit- 
ted curves are somewhat too shallow in the center of the 
minimum, though the effect is hardly larger than the ob- 
servational error margins. A physical explanation may be 
found in the treatment of the reflection effect. Our Wilson- 
Devinney-based solution code uses bolometric albedos and 
black body radiation to characterize the reflectivity of 
the irradiated surface layers. The upper limit of albedo 
is therefore restricted to 100% for complete reradiation 
of the incident flux, and in our case the adjusted albedo 
of the cool secondary was actually limited by this con- 
straint. Albedos larger than 1 would be "unphysical" , if 
local non-Planckian deviations of the intensity distribu- 
tion are not accounted for. Removing the constraint on 



A2 would result in an increased adjusted value of A2 of 
^ 1.1, and a deepening of the secondary minimum by a 
few millimag. Because of the tininess of this effect and 
due to the highly complex physical problems connected 
with the external irradiation of stellar atmospheres we pre- 
ferred not to use such purely phenomenological approxi- 
mation and remained with the classical reflection effect 
approach. Also, a suspected shift of the secondary mini- 
mum in the B curve by a few 10^^ phase units relative 
to phase 0.50 cannot be modeled without introducing ad- 
ditional free parameters. If real, it would conceivably be 
caused by a non-axisymmetrical brightness distribution of 
the heated surface area of the secondary component. 

The 1 a standard deviation of the normal points from 
the computed light curve amounts to about 0.004 mag. It 
should be noted that the best fit in the B passband (Agg 
= 435 nm) requires a third light contribution of 3.7%, 
while the R curve (Aog — 641 nm) is fully consistent with 
= 0. The same effect was observed in the closely re- 



lated sdB binary HW Vir (Wood et al. |1993D , in which 
the U and B filter measurements deviate from the com- 
puted curves resulting from a simultaneous UBVR fit by 
a similar amount, whereas the R curve shows no such sys- 
tematic effect. A possible explanation might arise from the 
strong reflection effect in both systems: though energy is 
conserved in the bolometric reradiation of the heated sur- 
face layers, there could be some frequency redistribution 
of photons inherent in the complex scattering process in 
a way that higher-energy photons have larger scattering 
cross sections than those incident with longer wavelengths 
and/or are transformed into lower-energy photons by var- 
ious absorption and reemission processes. 

The photometric solution yields a close detached sys- 
tem configuration for HS 0705+6700. A meridional inter- 
section of the surface equipotentials and the inner critical 
Roche lobe is shown in Fig. |^. The cool secondary is found 
to be relatively closer to its Roche lobe than the primary, 
which appears only slightly distorted by tidal interaction 
(''point /?'poic ~ 1-03 and 1.09, respectively, for the primary 
and secondary components). 3D-simulated aspects of the 
system at various orbital phases under a viewing angle 
corresponding to the derived orbital inclination of 84?4 
are displayed in Fig. ^. 

Though the simplex algorithm is a more powerful and 
numerically stable parameter optimization method than 
the differential corrections procedure, which is usually 
used in combination with the Wilson-Devinney approach, 
it also cannot completely rule out the possibility of con- 
vergence into a local minimum in the multi-dimensional 
parameter space. For this reason we experimented with a 
multitude of different start parameter sets to verify that 
our final solution corresponds to the deepest global min- 
imum. We found good evidence that this solution indeed 
turned out to be preferred in most of our trial runs. We 
should, however, mention that there were a few other solu- 
tions which differ from our best fit mainly in their values 




Fig. 3. Aspects of HS 0705+6700 at orbital phases 0.25, 0.45, 0.50, 0.575 (left column, from top to bottom), 0.75, 
0.95, 0.00, 0.075 (right column), viewed under the orbital inclination angle of 84?4. 



of the mass ratio q. While the best fit gave q — 0.278, 
other solutions with larger alternative values oiq — 0.346, 
0.460, and 0.594 were also found. Yet the standard devia- 
tions of the fits for q = 0.46 and 0.59 are worse by about 
5-10%, and they are also not consistent with our spec- 
troscopic results (inconsistent value of the mass function, 
see Sect. The photometric solution for q = 0.35 is of 
nearly identical numerical quality like our preferred one 
for q = 0.28, but in order to match the observed value 
of the mass function of f{m) — 0.0063 M©, the primary 
mass had to be assumed as < 0.3 M©. Hence a value of 
M2 < 0.1 M© would result. According to low-mass star 



models (e.g., Dorman et al. 1989), the mass-radius rela- 



tion would then yield a secondary radius of the order of 
0.10 — O.12i?0, which would be incompatible with the pho- 
tometrically derived value. 



It is a common problem inherent to light curve solu- 
tion techniques that it is not easy to identify the deep- 
est global minimum in the multi-dimensional parameter 
space. Since several parameters used to model the light 
curves are highly correlated, often a multitude of nearly 
equally good solutions can be found for different solution 
parameter sets. In our case, we also end up with vari- 
ous solutions with a comparable goodness of fit, which 
mainly differ in their q values (see above) . The light curve 
of HS 0705-^6700 is very simila r to those of HW Vir 
and PG 1336-018. Wood et al. (|l993|) and Kilkenny et 
al. ( 1998 ) both encountered the same problem when solv- 
ing the light curves of these systems and were able to 
obtain good solutions for a broad range of mass ratios. 
The simplex parameter optimization algorithm incorpo- 
rated in our MORO solution code is already a powerful 
tool with respect to the ability to automatically scan the 
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Table 3. Light curve solution of HS 0705+6700 



Fixed parameters: 



a;i(B) 
a:i(R) 
52-^ 



1.0 

1.0 

1.0 

0.32 

0.26 

0.19 

0.0 



Adjusted parameters: 



g(=M2/Mi) 

^2 

Teff(l) 
Teff(2) 

Li(B) ^ 

^3(B) ' 

Li(R) 

^3(R) ' 

X d 
Ol 



84?4 + 0?3 
0.278 + 0.019 
3.764 + 0.066 
2.606 + 0.051 
29 600K+800K 
2 900 K+600 K 
0.99997 + 0.00003 
3.7% + 0.6% 

0.99969 + 0.00031 
0.2 % + 0.2 % 

0.42 + 0.12 
0.62 + 0.17 
0.025 + 0.011 



Roche radii: ^ 



ri(pole) 


0.279 + 0.004 


ri (point) 


0.287 + 0.002 


ri(side) 


0.283 + 0.003 


ri (back) 


0.285 ± 0.003 


r2 (pole) 


0.219 + 0.001 


r 2 (point) 


0.238 + 0.004 


r2 (side) 


0.225 ± 0.002 


r2 (back) 


0.239 + 0.005 



^ bolometric albedo 

^ gravitational darkening exponent 

linear limb darkening coeffi cient; theoretical value taken 
from Wade & Rucinski(|l98|) 



199! 



radiation pressure parameter, see Drechsel et al. 
" relative luminosity Li/{Li +I/2); L2 is not independently 

adjusted, but recomputed from r2 and Teg (2) 
* fraction of third light at maximum 

^ fractional Roche radii in units of separation of mass cen- 
ters 



multi-dimensional space of adjustable parameters for the 
deepest minimum within broad parameter ranges. How- 
ever, the dependency of the goodness of fit on the chosen 
parameter set can be highly complex, or the relation is lo- 
cally "flat-bottomed" , so that a unique solution does not 
exist. 

To demonstrate this behaviour a number of solutions 
with different mass ratios were enforced by fixing q at 
various values in the surroundings of our best solution. 



Table 4. Goodness of light curve fit for various fixed q 
values compared with a of best fit solution with adjusted 

q 



q 




CTR 


O"tot 


0.20 


0.00433 


0.00429 


0.00428 


0.25 


0.00426 


0.00426 


0.00423 


0.278 1 


0.00407 


0.00403 


0.00401 


0.30 


0.00411 


0.00414 


0.00409 


0.35 


0.00408 


0.00402 


0.00401 


0.40 


0.00421 


0.00398 


0.00406 



adjusted q of final best solution 



The goodness of fit is measured by the standard deviation 
CTfit, which is defined by 

"fit — v^n J 

n - TO l^i^i Wi 
where n is the number of observations (normal points) 
and TO the number of adjustable quantities; O,; — Ci are 
the residuals between normal points and fitted curve, and 
Wi are the individual weights of the normal points. The 
simplex routine uses these same am values also internally 
during each iteration step to control the optimization pro- 
cedure. Table ^ summarizes the final atit values of conver- 
gent solutions obtained for a grid of arbitrary q values, 
which were kept fixed during these trial runs, while the 
rest of light curve parameters was adjusted to yield the 
best possible numerical representation of the observed B 
and R curves. We give tTfit values for the composite set of 
simultaneously adjusted B and R light curves, (Ttot, as well 
as the individual am values of the B and R curves, ctb and 
(Tr. The results for the trial solutions at q = 0.20, 0.25, 
0.30, 0.35, and 0.40 are compared with our selected best 
solution for q = 0.278, which was obtained by adjusting q 
together with the whole parameter set. It is evident that 
there are two ranges of q around 0.28 and 0.35, where the 
standard deviations occupy local minima in the parameter 
space. The trial solution for q = 0.35 (fixed) nearly coin- 
cides with one of our previously discussed unrestricted so- 
lutions with q treated as a free parameter (resulting in q — 
0.346). The distinction between the two alternative q val- 
ues 0.28 and 0.35 apparently cannot be made solely based 
on numerical standards, since the difference of the respec- 
tive a values is too small. Fortunately, in our case the 
available spectroscopic information (mass function, mass- 
radius relation for the low-mass dwarf companion) could 
be used to sort out inconsistent photometric solutions. 
The small error given in Table | for g = 0.278 ± 0.019 
simply reflects the formal fit uncertainty within the lo- 
cal minimum in parameter space around q — 0.28, and 
does not account for systematic effects like ambiguity of 
solutions. 
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4. Spectroscopic analysis 



Table 5. Atmospheric parameters of HS 0705+6700 



The optical spectra allow to measure the radial velocity 
curve with good phase coverage as well as a quantitative 
spectral analysis to determine the atmospheric parameters 
and the projected rotational velocity. 

Radial velocity curve 

Radial velocities from all Twin spectra were derived by 
cross-correlating the observed helium and Balmer line 
spectra to a synthetic spectrum calculated from a model 
atmosphere (see below) and are listed in Table ^ The 
measurements are accurate to about ±15kms^^ (Sc). 

Since this system is single-lined the analysis is straight- 
forward. A comparison of the best fit radial velocity curve 
and the measured values is shown in Fig. ^. The semi- 
amplitude is Ki = 85.8 ± 3.7kms^^and the systemic ve- 
locity 70 = -36.4 ±2.9 km s^i. 

From the period and semi-amplitude the mass function 
is derived: /(m) = 0.00626 ± 0.00081 Af©. 

4-2. Spectrum fitting 

In order to improve the S /N ratio the Twin spectra were 
radial velocity-corrected and then coadded. The coadded 
as well as the individual spectra were analyzed to derive 
atmospheric parameters using NLTE model atmospheres 
(Napiwotzki 1997 ) and a procedure described by Napi- 
wotzki et al. ( 119991) . 

Matching the synthetic Balmer (H/3 to He) and He I 
(4026 A, 4471 A and 4922 A) line profiles to the obser- 
vations resulted in a determination of effective tempera- 
ture, gravity and He abundance. The resulting fit for the 
coadded spectrum is displayed in Fig. || and gave Toh ~ 
28 800 ± 300 K, logg = 5.40 ± 0.04 and log(nHc/riH) = 
—2.68 ± 0.05. The quoted errors are derived from the 
procedure applied to fit the line profiles. 

We also analyzed the individual spectra in the same 
way as described above. Two spectra (at phases 0.701 and 
0.855) were not useable for a quantitative spectral analy- 
sis because of insufficient S/N. These two phases, however, 
are covered by better spectra (at phases 0.711 and 0.849) 
anyway. The results for the remaining 13 phases are sum- 
marized in Table ||, and are plotted in Fig. ^. Because of 
the reflection effect the Balmer lines might be distorted by 
reflected light from the secondary. Indeed, in the case of 
HW Vir Wood & Saffer ( |1999D found that Tcft varied by 
1500 K and log 5 by 0.1 dex around the orbit. Two spec- 
tra of HS 0705+6700 were taken near primary and another 
near secondary minimum, for which we do not expect any 
contamination by reflected light. 

The solution of the B light curve indicated the pres- 
ence of third light. This could cause a systematic error 
for the spectroscopic analysis. The spectral characteristics 
of this light, however, are unknown. We made an exper- 



phase 


Teff (K) 


log 3 




-log(He/H) 


0.107 


27300 ± 836 


5.41 ± 


10 


2.50 + 0.11 


0.203 


27720 ± 858 


5.29 ± 


11 


2.82 + 0.14 


0.276 


29748 ± 901 


5.48 ± 


14 


2.60 + 0.22 


0.303 


29063 ± 697 


5.30 ± 


10 


2.44 + 0.13 




98079 + R'\A 


5.28 + 


11 


9 + n 1 9 


0.428 


28947 ± 1600 


5.12 + 


23 


2.99 + 0.49 


0.505 


28414 ± 1002 


5.30 + 


13 


2.63 + 0.12 


0.574 


29928 ± 828 


5.43 + 


14 


2.81+0.23 


0.601 


28649 ± 751 


5.59 + 


11 


2.54 + 0.11 


0.711 


30126 ± 649 


5.51 + 


11 


2.89 + 0.19 


0.742 


28813 ± 849 


5.50 + 


12 


2.78 + 0.15 


0.849 


28186 ± 842 


5.26 + 


10 


2.58 + 0.11 


0.992 


27102 ± 1009 


5.37 + 


13 


2.94 ± 0.20 


Mean: 


28621 ± 919 


5.37 + 0.12 


2.70 + 0.17 


coadded 


28755 ± 280 


5.40 + 0.04 


2.68 ± 0.05 



iment on an individual spectrum and subtracted 3% of 
the continuum from the blue spectra and repeated the fit. 
In this case the effective temperature is lowered by about 
the same amount as the la errors listed in Table |. The 
change in gravity is only a few hundredth of a dex, less 
than the errors listed in this table. As possibly indicated 
by Fig. ^, there might be a slight variation of Tes with 
orbital phase, but not so for loggr. Our spectra are prob- 
ably of insufficient quality to reveal any phase-dependent 
variations of the atmospheric parameters at such low level. 

We finally adopted Teff = 28 800 ± 900 K, log 5 = 
5.40 ± 0.1 and log(nHc/?^H) = -2.68 ± 0.15 for the at- 
mospheric parameters of HS 0705+6700. The given errors 
were estimated from the scatter of the Teff and logg val- 
ues derived from the individual spectra and hence would 
cover the bandwidth of any possibly present systematic 
variation with phase. 

We finally note that the spectroscopic Toff is in excel- 
lent agreement with the one following from the photomet- 
ric solution (Sect. |3.2| ). 

4.3. Projected rotational velocity 

The line profiles of Ha and He 1 6678A could be used to 
determine the projected rotational velocity. We calculated 
the Ha profiles from the flnal model and convoluted it 
with rotational profiles for various rotational broadening 
parameters. The best fit was again obtained by a proce- 



dure as described by Heber et al. (1997), which resulted in 
a projected rotational velocity of usinz = UOti^kms-i 
(3(7 errors). The corresponding fit is displayed in Fig. ^ 
The rotation of HS 0705+6700 is expected to be tidally 
bound. In this case we can estimate the projected rota- 
tional velocity from the orbital period and the radius (see 
Sect. H 

result from spectrum synthesis 



to be 12lli[3kms ^ in good agreement with the 




Fig. 4. Radial velocity curve of the sdB primary component in a circular orbit; velocities are heliocentrically corrected; 
the semi-amplitude Ki is 85.8 km s~^, the systemic velocity 70 = — 36.4kms~^; error bars correspond to 3 a margins. 



HS 0705 + 6700 






AA/A 



Fig. 5. Fit of the Balmcr and helium lines in the coadded 
spectrum of HS 0705+6700 by synthetic spectra calculated 
from NLTE model atmospheres. 



5. System parameters 

5.1. HS 0705+6700 

Even though the radial velocity amplitude of the compan- 
ion is unknown we have sufficient information to calculate 
absolute system parameters due to the M — R relation 
provided by the logg determination (Sect. |4.2| ). 

From the light curve we have derived a period of 
8263.87 s, an inclination angle of i = 84?4, and a mass 
ratio q — 0.278. Furthermore, we determined the radii for 
the binary components Ri/a = 0.284 and R2/a = 0.230, 
where a is the separation of their mass centers. Using our 
measured Ki — 85.8 km s^^, log 171 = 5.40, and the in- 
clination i, we derive Mi = 0.483 Mq, M2 = 0.134Mq, 
Ri = 0.230 i?o, R2 = 0.186 i?0, and a = 5.65 • 10^° cm 
(= 0.81 i?0). From these numbers we conclude that the 
companion is an M dwarf. According to the models for 



low-mass dwarf stars of Dorman et al. (1989), a star of 
0.14 M© would have a radius of 0.17 i?©, fully consistent 
with the observations. 

The above given parameters correspond to the photo- 
metric solution for q = 0.278. As can be seen from Fig. ||, 
the position of HS 0705-1-6700 in the HRD is well compati- 
ble with those of other EHB stars as well as with EHB evo- 
lutionary tracks of stars with corresponding masses. Con- 
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0.4 0,6 
Orbital phase 



Fig. 6. Variation of atmospheric parameters with phase; 
the error bars result from the fit procedure described 
by Napiwotzki et al. 



1999 



sidering the alternate photometric solution for q = 0.35, 
one would have to assume a lower primary mass of say 
Ml = 0.3 M©. In principal, such mass could result as end 
product of the common envelope evolution. The HRD po- 
sition of the primary might then also match a track leading 
to a helium white dwarf. However, as already mentioned 
earlier, the reduced mass of the M dwarf would require 
a smaller radius for this star, hardly compatible with the 
light curve solution. 



Fig. 7. Determination of the projected rotational velocity; 
synthetic Ha and He I line profiles calculated from the fi- 
nally adopted model and the best fit v sini = 110kms~^ 
are compared to the coadded spectrum; the same line pro- 
files predicted for a non-rotating star are also shown as 
dashed lines. 



Given the most probable parameters and system con- 
figuration of HS 0705-f6700 (Fig. |), the further evolu- 
tion might lead to the formation of a cataclysmic system 
with a period below the period gap. The process of an- 
gular momentum loss by gravitational radiation would re- 
quire about 10^ years for the binary components to merge. 
Should the radius of the sdB star along its track towards 
the white dwarf cooling sequence, however, increase suffi- 
ciently to reach its limiting lobe, another episode of rapid 
mass transfer and loss of mass and angular momentum, 
perhaps even another common envelope phase could oc- 
cur within a much shorter time scale. Yet, such prospects 
should be subject of detailed evolutionary calculations ly- 
ing beyond the scope of this paper. 

5.2. Similar systems 

HS 0705+6700 is only the third eclipsing binary known to 
consist of an sdB star and an M dwarf companion. An or- 
bital period of about 0.096 days places HS 0705+6700 in a 
period range, which coincides with the period gap of cata- 
clysmic variables. There are two other surprisingly similar 
systems - HW Vir (P f« 0.117 d) and PG 1336-018 (P « 
0.101 d, ) - which have equally short periods within the 
CV p eriod gap between about 2 and 3 hours (e.g., Shaffer 



1992). A non-eclipsing system (PG 1017-086) with an even 
shorter period (P = 0.073 d) than HS 0705+6700 was dis- 
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Table 6. System parameters of three similar eclipsing sdB binaries 





xi vv V ir 


rVj looD-Ulo 


xlo UtUo-tOtUU 


Author 


Wood et al. 1993 


Kilkenny et al. 1998 


this paper 


data 


UBVR 


UVR 


BR 


P 


2h 


2*^ 25"" 




q = M2/M1 


~ 0.3 


0.3 


0.278 


i 


80°. 6 


81° 


84°. 4 


Ki 


87.9 km s"^ 


78 kms"^ 


85.8 km s~^ 


Ti 


- 33 000 K 


33 000 K 


29 600 K 


T2 


~ 3 700 K 


~ 3 000 K 


2 900 K 


log 51 


5.64 


5.7 


5.40 


Ml 


0.54 Mq 


- 0.50 Mq 


0.483 Mq 


M2 


0.18 M© 


0.17 Mq 


0.134 Mq 


Ri 


0.183 Rq 


0.165 i?0 


0.230 Rq 


R2 


0.188 Rq 


0.175 _Rq 


0.186 Rq 


a 


0.89 Rq 


0.79 Ti© 


0.81 i?Q 



0.471 V,, « 




1: HW Vir 

3: PG 1336-018 

3: PG 1017-086 



30 

T,„/1000K 



Fig. 8. Comparison of the position of HS 0705+6700 in 
the (Teff, logg) diagram to that of HW Vir, PG 1336- 
018 and PG 1017-086, EHB stars studied by Maxted et 



al. ( 2001a ) and to evolutionary EHB tracks of Dorman et 
al. ( IT993I) ; the tracks are labeled with the stellar mass. 



covered recently by Maxted et al. ( pOOlb ). The shrinkage 
of the orbit from an originally much wider binary with 
a red giant component to a component separation of less 
than one solar radius can be explained by a common enve- 
lope stage, during which a large fraction of orbital angular 
momentum is lost. The system parameters of the three re- 
lated sdB binaries are compared in Table ^. 

In Fig. ^ we compare the position of the four systems 
with a sample of sdB stars recently analyzed by Maxted 



et al. ( 2001a ) and with calculations for EHB evolution. 
As can be seen all four stars lie in the same region as the 
other sdB stars on the EHB band, and are therefore core 
helium burning stars with very small inert hydrogen-rich 
envelopes. In order to evolve to their present configuration 
the systems must have gone through a common envelope 
phase, when the primary was near the tip of the red gi- 
ant branch. Both the very short periods as well as the 
extremely thin hydrogen layers on top of the helium cores 
are conceivably explained by such common envelope evo- 
lution connected with extensive mass loss. 

In the case of HS 0705+6700 and the mentioned related 
systems the evolutionary scenario through which the sdB 
stars were formed is most probably the same. This raises 
again the question whether all sdB stars are members of 
close binary systems and originate from the same common 
envelope evolution. To improve the statistical significance 
of such a hypothesis, it will therefore be of major impor- 
tance to search for and analyze a larger number of sdB 
candidates with special emphasis of their possible bina- 
rity 



6. Conclusions 

The sdB star HS 0705+6700 was discovered to be an 
eclipsing binary with an orbital period of 8263.87 s. The 
analysis of the light curve revealed that the mass ra- 
tio of the system is g = 0.28 and the inclination 84?4. 
The companion does not contribute to the optical light of 
the system, except for the light reflected from the hemi- 
sphere facing the sdB star. The semi-amplitude of the ra- 
dial velocity curve Ki — 85.8 km s^^ and a mass function 
of /(to) = 0.00626 Mq were derived. Accurate absolute 
masses and radii of the sdB primary and the M dwarf 
companion could be determined. The spectroscopic analy- 
sis of the sdB component resulted in Toff = 28 800 ±900 K, 
logg = 5.40 ± 0.1, and log(nHo/nH) = -2.68 ± 0.15. 
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The position of the sdB star in the log g - Teft diagram 
coincides with the domain of many other EHB stars. The 
current location can be matched by the evolutionary track 
for such a star with corresponding mass. The photometric 
value of the mass ratio oiq = 0.28 and the respective mass 
of the secondary are in good agreement with the observed 
mass function and with the predicted mass-radius relation 
for low-mass dwarf stars. Based on these arguments an al- 
ternate photometric q value of 0.35 appears less probable. 

The further evolution of the system might lead to the 
formation of a cataclysmic system with a period below the 
CV period gap. 
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